arXiv:nucl-th/9706083vl 28 Jun 1997 


Importance of Non-Linear Couplings in Fusion Barrier 
Distributions and Mean Angular Momenta 

K. Haginof N. Takigawa^, M. Dasgupta^, D.J. Hinde^, and J.R. LeigD 

^ Department of Physics, Tohoku University, Sendai 980-77, Japan 

^ Department of Nuclear Physics, 

Research School of Physical Sciences and Engineering, 

Australian National University, Canberra, ACT 0200, Australia 

Abstract 

The effects of higher order coupling of surface vibrations to the relative motion on 
heavy-ion fusion reactions at near-barrier energies are investigated. The coupled channels 
equations are solved to all orders, and also in the linear and the quadratic coupling 
approximations. It is shown that the shape of fusion barrier distributions and the energy 
dependence of the average angular momentum of the compound nucleus can significantly 
change when the higher order couplings are included. The role of octupole vibrational 
excitation of in the -I- fusion reaction is also discussed using the all order 

coupled-channels equations. 

1. Introduction 

The analysis of the fusion process in terms of the barrier distribution has generated 
renewed interest in heavy-ion collisions at energies below and near the Coulomb barrier]^. 
In a simple eigenchannel approach, couplings of the relative motion of the colliding nuclei 
to nuclear intrinsic degrees of freedom result in the single fusion barrier being replaced by 
a distribution of potential barriers. A method of extracting barrier distributions directly 
from fusion excitation functions was proposed[0, and stimulated precise measurements 
of the fusion cross sections for several systems. The analysis of the barrier distribution 
clearly shows the effects of couplings to several nuclear intrinsic motions in a way much 
more apparent than in the fusion excitation function itself]^]. 

Theoretically the standard way to address the effects of the coupling between the rela¬ 
tive motion and the intrinsic degrees of freedom is to solve the coupled-channels equations, 
including all the relevant channels. Most of the coupled channels calculations performed 
so far use the linear coupling approximation, where the coupling potential is expanded in 
powers of the deformation parameter, keeping only the linear term. Whilst this approach 
reproduces the experimental data of fusion cross sections for very asymmetric systems. 
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it does not explain the data for heavier and nearly symmetric systems^]. Thus, it is 
of interest to examine the validity of one of the main approximations in these calcula¬ 
tions, namely the linear coupling approximation, and see whether the effects of non-linear 
coupling improve the agreement between data and the theoretical calculations for such 
systems. Even in asymmetric systems, the non-linear couplings might be important to 
reproduce precisely measured data. 

In this contribution, we report the results of the coupled channels calculations which 
include the couplings to all orders]^, |^. The results of these calculations for fusion cross 
sections, average angular momenta and barrier distributions are compared with those 
using the linear and the quadratic coupling approximations. It is seen that the linear 
coupling approximation is not valid even in systems where the coupling is weak, and that 
higher order couplings strongly influence the calculated barrier distributions. 

2. Coupled channels equations in all orders 


Consider the problem where the relative motion between colliding nuclei couples to 
a vibrational mode of excitation of the target nucleus. For simplicity excitations of the 
projectile are not considered in this section. It is straightforward to extend the formulae 
to the case where many different vibrational modes are present and where projectile 
excitations also occur. For heavy ion fusion reactions, to a good approximation one can 
replace the angular momentum of the relative motion in each channel by the total angular 
momentum J[§|. This approximation, often referred to as no-Coriolis approximation, will 
be used throughout this paper. The coupled channels equations then read 

^nir) + ^nmir)^/jmir) = 0 , ( 1 ) 

m 

where r is the radial component of the coordinate of the relative motion, /i the reduced 
mass, Ecm the bombarding energy in the center of mass frame and hu is the excitation 
energy of the vibrational phonon. 14im are the coupling matrix elements, which in the 
collective model consist of Coulomb and nuclear components. 

For the Coulomb coupling, it has been shown that inclusion of up to the hrst order 
term is sufficient 0. The matrix elements of the Coulomb coupling in Eq. (1) thus read 
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( 2 ) 


where is the amplitude of the zero point motion. It is related to the deformation 
parameter Px by ao = + 1|01- 

In the collective model, the nuclear interaction is assumed to be a function of the sep¬ 
aration distance between the vibrating surfaces of the colliding nuclei. It is conventionally 
taken as 


V^^\r,axo) 


_H_ 

1 -F exp[(r - Rp- Rt - \f^^RTa\o)/a] 


( 3 ) 


where the surface coordinates are related to the phonon creation and annihilation 
operators by cta^ = ao{axfj, + {—)^ax-fj^). Volume conservation introduces a small term 
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which is non-linear with respect to the deformation parameter axo in the denominator of 
the above Eq.(3). This is ignored for simplicity in the present study. 

The conventional nuclear coupling form factor in the linear coupling approximation is 
obtained by expanding Eq.(3) with respect to axo and keeping only the linear term. It is 
given by 


= VN{r)6n,ni “ iy/n6n,m+i + Vn + 16 n,ni-l), 


( 4 ) 


where VAr(r) = — Vo/[l -|- exp((r — Rp — i?r)/a)] is the nuclear potential in the entrance 
channel. 

In order to take into account the effects of the couplings to all order, we evaluate 
the nuclear coupling matrix elements without introducing the expansion. Denoting the 
eigenvalue of axo by the matrix elements of the nuclear coupling then read 



dxU*^{x)Um{x) 


_ 

1 -|- exp[(r — Rp — Rp — 


\/^RTx)/a\ 


( 5 ) 


Here Un{x) is the eigen function of the n -th excited state of the harmonic oscillator and 
is given by 


Un{x) = 


H„ 


X 


”VV2«o 


^-x'^lAal 


where Hn{x) is the Hermite polinomial with rank n. 
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3. Validity of the linear coupling approximation 


3.1. Nearly Symmetric Systems 


We now present the results of our calculations of fusion cross sections, average angular 
momenta of the compound nucleus, and fusion barrier distributions. We hrst discuss 
heavy nearly symmetric systems. We analyse in particular ^'^Ni -|- ®®Zr reactions which 
are typical examples where the conventional coupled channels calculations with the linear 
coupling approximation fail to reproduce the fusion cross sections and average angular 
momentum data^]. Our aim is to investigate whether the failure is due to the linear 
coupling approximation by performing linear, quadratic and full coupling calculations. 

We take into account the couplings up to two phonon states of the quadrupole sur¬ 
face vibrations of ®‘^Ni, and of the octupole vibration of ®®Zr. We also take their mutual 
excitations into account. We ignore the effects of transfer reactions, because it has been 
reported in Ref. |Q] that they have only small effects on the fusion cross sections and the 
average angular momenta in this reaction. The excitation function of the fusion cross 
section for this reaction obtained by numerically solving the coupled channels equations 
is compared with the experimental data in Fig. 1 (upper panel). The experimental data, 
taken from Ref. |^, consist only of the evaporation residue cross sections, and do not 
include hssion following fusion. The dotted line is the results in the one dimensional 
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potential model, i.e. without the effects of channel coupling. The dot-dashed line is the 
results of the coupled channels calculations when the linear coupling approximation is 
used. They considerably underestimate the fusion cross sections at sub-barrier energies. 
The situation is slightly improved when the quadratic coupling approximation^ is used, 
i.e. when the nuclear coupling potential up to the second order of the deformation pa¬ 
rameter, is included (dashed line). However, there still remain considerable discrepancies 
between the experimental data and the results of the coupled channels calculations. When 
we include couplings to all order, we get the solid line, which agree very well with the 
experimental data. Dramatic effects of the higher order couplings on fusion cross sections 
are observed, especially at low energies. 

The lower panel in Fig. 1 compares the results of our calculations of the average 
angular momentum of the compound nucleus with the experimental data as a function 
of the bombarding energy. We again observe that the experimental data are much better 
reproduced by taking the effects of couplings to all orders into account. We thus conclude 
that coupling to all orders are essential to simultaneously reproduce the fusion cross 
sections and the average angular momentum data for heavy (nearly) symmetric systems. 
This is in agreement with the calculations required to fit the barrier distribution for ^®Ni 
- 1 - react ion ||l0|. 

3.2. Very Asymmetric Systems 

We next consider the effects of higher order couplings for very asymmetric systems 
where the product of the charges ZpZx is relatively small. For such systems, the coupled 
channels calculations in the linear coupling approximation have achieved reasonable suc¬ 
cess in reproducing fusion excitation functions. However, no study has been performed to 
see whether the effects of higher order couplings on the angular momentum distribution 
of the compound nucleus and on the barrier distributions are small. In this subsection 
we re-analyse the experimental data for the -|- ^^^Sm reaction, for which the effects 
of couplings to phonon states on fusion barrier distributions were shown experimentally 
for the first time[]^. 

The fusion barrier distribution for this system is intimately related to the octupole 
vibration of ^^^Sm, and that the quadrupole vibration plays only a minor role[|^. Ac¬ 
cordingly, we ignore the effects of the couplings to the quadrupole phonon states of ^^^Sm 
and include only the single octupole phonon state at 1.81 MeV. For simplicity in the 
calculations we ignore excitation of the projectile. These effects will be discussed in Sec. 
4. 

The upper panel of Fig. 2 shows the experimental fusion excitation function taken 
from Ref. |l|, and the theoretical calculations. The meaning of each line is the same as 
in Fig.l. We observe that the agreement of the theory and experiment appears to be 
improved only slightly by the inclusion of coupling to all orders. The barrier distribution, 
on the other hand, reveals significant changes due to the higher order couplings (see the 
lower panel of Fig. 2). Comparing the results of the linear coupling approximation (the 
dot-dashed line) with those of the all order coupling (the solid line), one observes that the 
higher order couplings transfer some strength from the lower barrier to the higher barrier, 
and at the same time lower the peak position of both barriers. 
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Figure 3 shows the coupling matrix element between the ground state and the one 
phonon state Voi(r) as a function of the inter nuclear separation distance r. The dotted 
line is the coupling matrix element in the linear coupling approximation, while the solid 
line includes the coupling to all orders. One can see that the linear coupling approximation 
underestimates the coupling strength at the barrier position of the uncoupled barrier 
around r= 10.8 fm. On the other hand, it overestimates the coupling strength in the 
inner region around r=8.5 fm. As we will see in the next section, the latter fact is 
important in discussing the effects of the projectile excitation. 

The effects of higher order couplings become more signihcant when there exist more 
than two channels. In order to demonstrate this, we show in £g. 4 the barrier distributions 
for the ^®0 + reaction, where the double octupole phonon excitations are allowed 

in the harmonic limit. One observes dramatic effects of higher order couplings on the 
fusion barrier distribution. It is thus clear that high precision measurements should be 
analysed using all order couplings even when coupling is weak as a result of the small 
charge product. The detailed studies on the effects of double phonon couplings including 
anharmonic effects in this reaction is reported separately BH- 

4. Role of projectile excitation in fnsion 


Let us now discuss the role of the octupole vibrational excitation of in the + 
fusion reactions[^. Contradictory conclusions have been reported regarding the 
role of projectile excitation in the fusion reactions between and samarium isotopes. 
Calculations of the fusion cross-section for the -|- ^^^Sm reactions in ref. [|1^ indicated 
the importance of the excitation of In marked contrast, no specihc features appear in 
the measured barrier distribution for the reactions which can be associated 

with the excitation of rather, it was argued in Ref. that a good theoretical repre¬ 
sentation of the experimental fusion barrier distribution is destroyed when the projectile 
excitation is included. 

Both of these conclusions are based on the comparison of the experimental data with 
the results of simplified coupled-channels calculations, where the linear coupling approx¬ 
imation is used to describe the vibrational excitation of the projectile. In the previous 
section, we have shown that the linear coupling approximation is not valid even in systems 
with weak coupling, and that higher order couplings strongly influence the fusion barrier 
distribution. The coupling to the octupole vibrational state of is strong because of the 
large deformation parameter indicated by the strong E3 transition. It is therefore very 
likely that the fusion barrier distribution calculated in the simplified coupled-channel 
codes using the linear coupling approximation does not represent the true fusion barrier 
distribution. 

In order to test the validity of the simplihed coupled-channels calculations, here we 
hrst calculate the excitation function of the fusion cross section and the fusion barrier 
distribution in the linear coupling approximation. The results are shown in Fig. 5. The 
dotted line is the results when is treated to be inert. This calculation well reproduces 
the features of the experimental barrier distribution. The results of calculations including 
the excitation of the lowest-lying octuple state of are shown by the solid line. Though 
the experimental barrier distribution around the lower energy peak at (~ 60 MeV) is 
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reproduced, significant strength is missing around the higher energy peak near 65 MeV. 
A similar discrepancy between theory and experimental data was encountered in Ref. |jl|] , 
where calculations were performed using a modified version of the CCFUS code (the long- 
dashed line). Clearly both calculations which treat the octupole excitation of in the 
linear coupling approximation fail to reproduce the experimental barrier distribution. 

The results of coupled-channels calculations, where the couplings to the octupole vi¬ 
brations of both and are treated to all orders, are shown in figure 6. It is 

remarkable that these calculations re-establish the double-peaked structure seen in the 
experimental data, which was missing in the linear coupling calculations. Indeed, the 
barrier distribution obtained by including the coupling to the octupole vibration of 
to all order looks very similar to that obtained by totally ignoring it, apart from a shift 
in energy. A shift of 2 MeV (dashed line) of the former is required for the two calcu¬ 
lated distributions to coincide. This shift is consistent with the general conclusion that 
the main effect of the coupling to an inelastic channel whose excitation energy is larger 
than the curvature of the bare fusion barrier, i.e. an adiabatic coupling, is to introduce a 
static potential shift |^, and hence, the shape of the barrier distribution does not change 
unless the coupling is very strong and the coupling form factor itself has a strong radial 
dependence. The effects of these excitations can then be included in the ‘bare’ potential 
in the coupled-channels calculations. Thus for -|- ^^^Sm, where potential parameters 
for the calculations are obtained from a fit to the high energy data, the effects of oc¬ 
tupole vibration of are already included. The explicit inclusion of the coupling to 
octupole vibration then leads to double counting which manifests itself by introducing an 
additional shift in the barrier (or barrier distribution) as observed earlier. 


5. Summary 


We have shown that higher order couplings to nuclear surface vibrations play an 
important role in heavy ion fusion reactions. The inclusion of the coupling to all orders 
in coupled-channels calculations was shown to be crucial to reproduce the experimental 
fusion cross sections and the average angular momenta for heavy symmetric systems. 
We performed coupled-channels calculations also for the reactions as an 

example of very asymmetric systems where the coupling is weaker. It was found that 
higher order couplings to the vibrational states of the target nucleus result in a non- 
negligible enhancement of the fusion cross sections and a significant modification of barrier 
distributions even in very asymmetric systems. Our studies warn that spurious conclusions 
could be reached regarding the nature of couplings if high quality experimental data are 
compared with simplified calculations in the first order approximation. We also studied 
the role of projectle excitation in the fusion reactions. Using the calculations 

with full order coupling we have shown that the major effect of the excitation of the 
octupole state at 6.1 MeV of is to renormalize the static potential barrier without 
significantly modifying the shape of the barrier distribution. 
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Figure Captions 

Fig.l lExcitation function of the fusion cross section (upper panel) and the average angular 
momentum of the compound nucleus (lower panel) for the ®^Ni -|- ®^Zr reactions. 

Fig.2 lExcitation function of the fusion cross section (upper panel) and the fusion barrier 
distribution (lower panel) for the -|- ^^"^Sm reactions. In the coupled channels calcula¬ 
tions, the projectile is assumed to be inert, while the single octupole phonon state of the 
target nucleus is taken into account. 

Fig.3 :The coupling matrix element between the ground state and the one phonon state for 
the reaction as a function of the separation distance between the projectile 

and target.. 

Fig.4 :Same as the lower panel of £g. 2, but for the case where the double octupole 
phonon excitations of are included in the harmonic limit. 

Fig.5 : Effects of the projectile excitation on the fusion reactions. The linear 

coupling approximation is used in the coupled-channels calculations. In all calculations, 
the effects of the octupole vibration of ^^^Sm are taken into account. 

Fig.6: Same as Fig.5, but for the case where the coupled-channels calculations have been 
performed to all order coupling. 
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